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ABSTRACT 

Models of galaxy formation predict that gas accretion from the cosmic web is a primary driver of 
star formation over cosmic history. Except in very dense environments where galaxy mergers are also 
important, model galaxies feed from cold streams of gas from the web that penetrate their dark matter 
haloes. Although these predictions are unambiguous, the observational support has been indirect so 
far. Here we report spectroscopic evidence for this process in extremely metal-poor galaxies (XMPs) 
of the local Universe, taking the form of localized starbursts associated with gas having low metallicity. 
Detailed abundance analyses based on Gran Telescopio Ganarias (GTG) optical spectra of ten XMPs 
show that the galaxy hosts have metallicities around 60% solar on average, while the large star-forming 
regions that dominate their integrated light have low metallicities of some 6% solar. Because gas mixes 
azimuthally in a rotation timescale (a few hundred Myr), the observed metallicity inhomogeneities 
are only possible if the metal-poor gas fell onto the disk recently. We analyze several possibilities for 
the origin of the metal-poor gas, favoring the metal-poor gas infall predicted by numerical models. 

If this interpretation is correct, XMPs trace the cosmic web gas in their surroundings, making them 
probes to examine its properties. 

Subject headings: galaxies: abundances - galaxies: dwarf - galaxies: evolution - galaxies: formation 
- galaxies: structure - intergalactic medium 


1. MOTIVATION 


Metals are primarily produced by stars. Thus, 
metal-poor objects are chemically unevolved, providing 
a gateway to sample physical conditions and processes 
characteristic of early phases in the Universe. XMP 
galaxies, defined to have less than one tenth of the 


metals in the solar composition (iKunth fc OstlinI 1200011 . 
are rare in the local Un iverse (< 0.2% of all galaxies; 
iMorales-Luis et aT]l2011[l . They tend to be dwarf galax¬ 
ies b ut, unlike other dwarfs (jKobulnickv fc SkillmanI 
119961) . their chemical c o mposit io n may not be uni¬ 
form (iPanaderos et al.l 120061: lizotov et al.l 
Levesque et al.l 120111 : iSanchez Almeida et al.l 


2009; 


2013, 


2014bf) . with the low metallicity only in regions of in¬ 


tense star formation. Such non-uniformity is particularly 
revealing because the timescale for azimuthal mixing 
in normal galaxies is short, of the order of a fraction 
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of the rotational period (e.g ., Ide Avillez fc Mac Lowl 
120021: lYang fc KruinholS 12012 ). It would imply that 
the metal-poor gas in XMPs was recently accreted 
from a nearly pristine cloud, providing long-sought 
evidence for cosm ic accretion that has been diff icult to 
observe directly (|Sanchez Almeida et al.l l2014al) . For 
example, gas flows have been observed on the periphery 
of distant galaxies using absorption lines in the spectra 
of background sources. However, the material is often 
not significantly lower in metals than the galaxy itself, 
suggesting gas recycled from previous outflows, and 
therefore dy namically detached from the large-scale 
cosmic web (iLehner et all 120131: iGhurchill et al.! I2013t 
iLiang fc Chenll2014ll . 

Here we put on firm observational bases the existence 
of chemical inhomogeneities in most XMPs, thus un¬ 
equivocally showing their star formation to be feeding 
from external metal-poor gas. Measuring chemical inho¬ 
mogeneities in XMPs is technically challenging, requiring 
spectra of high quality from faint regions. We used the 
10.4 m GTCUH (Sect. [2]) to measure oxygen abundance 
along the major axes of ten XMP galaxie s (Tabled]) with 
a rob ust variant of the direct method (iPerez-Monterol 
12014 see Sect. [3]). In nine out of ten cases, sharp metal¬ 
licity drops were found (Sect.|4]). The origin of the metal- 
poor gas is analyzed in Sect. [SJ leaving as the most con¬ 
vincing hypothesis the gas infall from the web predicted 
by numerical models. 

2. OBSERVATIONS 


http://www.gtc.iac.es/GTChome.php 
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The ten galaxies in Table [J represe nt 7 % of all known 
XMP s when the project began (|Morales-Luis et ^ 
[MID, and they were chosen to be representative of the 
XMP family. They were observed with the long-slit mode 
of the spectrograph OSIRIS at GTC. We used an in¬ 
termediate resolution grism which renders the spectral 
range needed for abundance analysis in a single expo¬ 
sure (from 3700 A to 7000 A) with a dispersion of 7.8 A 
per pixel, enough to spectrally resolve the required emis¬ 
sion lines. The width of the spectrograph slit was set 
to 1 arcsec, which matches the typical seeing during ob¬ 
servation. Integration times amount to two hours per 
target. The spectra were reduced and calibrated us¬ 
ing standard modules of the package IRAFQ, with the 
line fluxes computed from the calibrated spectra using 
custom-made IDL routines. 

3. DETERMINATION OF METALLICITIES, 

STELLAR MASSES, AND STAR-FORMATION 
RATES 

HII regions are photoionized by young stars. Their 
emission-line spectrum allows the determination of the 
oxygen abundance independently of the details of the 
ionization source and the geometry of the nebula 
(|Osterbrocklll974 IStasihska et al.ll20I^ . The spectrum 
contains a number of collisionally excited emission lines 
that empirically provide electron-density and electron- 
temperature. Given temperature and density, atomic 
properties yield the expected emission per ion, therefore, 
the total flux of a particular emission line directly gives 
the total number of ions in the nebula responsible for 
such emission. The case of oxygen is particularly fa¬ 
vorable since, in HII regions, oxygen is in the form of 
0+ and 0++, and both species produce strong emis¬ 
sion lines. The total oxygen abundance results from 
adding the contribution of 0+ and 0++. A drawback 
of this technique (the so-called direct method; DM) is 
the need to measure fluxes of weak lines often at the 
noise level. Unfortunately, this problem affects the faint 
tails of our XMPs (Figs. [T] and H] ). We circumvent th e 
difficulty using HII-CHI-mistrjO (I Perez- Monteroll2014f) . 
a procedure that compares the brightest optical emis¬ 
sion lines with predictions of a grid of photoioniza¬ 
tion models covering a wide range of possible physical 
conditions. The resulting metallicity is calculated as 
the x^-weighted mean of the metallicities in the mod¬ 
els, with the x^-weighted standard deviation providing 
the error. This method leads to an oxygen abundance 
which is consistent with the DM determination within 
0.07 dex, i.e., within t he intrinsic uncertainty of the DM 
(|Perez-Monterol[201411 . The agreement implies that HII- 
GHI-mistry inherits the good properties of DM, so that 
the geometry and other assumptions needed to build 
the grid of models do not influence the inferred oxygen 
abundance. Provided that the emission lines are pro¬ 
duced by photoionization, the agreement between HII- 
GHI-mistry and DM guarantees that irrespectively of the 
source of ionizing photons, our abundance determina¬ 
tions along the XMP major axes are reliable and can be 
inter-compared. Photoionization is a good approxima¬ 
tion even for the faintest parts of XMPs, as attested by 

http://iraf.noao.edu/ 

http: //www. iaa. es/~epm/HII-CHI-mistry. html 
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Fig. 1.— Distribution of oxygen along the major axis of a typical 
XMP, J1132+57. Top: SDSS image with the red bars indicating 
the position of the slit of the spectrograph during observation, and 
the arrow pointing out the sense of growing distances in the bottom 
panel. Bottom: Variation of oxygen abundance (red symbols with 
error bars) and surface SFR (blue solid line) along the major axis 
of the galaxy as measured from emission-line spectra. The drop 
in abundance associated with the peak SFR coincides with the 
brightest knot of the galaxy. North and east are indicated by white 
arrows, and the bar on the top-left corner corresponds to 10 arcsec 
on the sky. 

the weakness of the line [NIIJA6583 A. In the cases where 
we have been able to measure it, its flux is less than a 
tenth of the Ha flux which, accordi ng to the BPT diag¬ 
nostic diagram (jBaldwin et al.llIQSip. discards ionizatio n 
by Active Galacti c Nuclei, shocks (iKewlev et al.ll20I3fl . 
and evolved stars l|Cid Fernandes et al.ll20lTfr 

The stellar masses of the galaxies given in Table [1] 
were determined by the MPA -JHU collaboration from 
galaxy-integrated magnitudes (|Brinchmann et ^ 12004 
ISalim et al.l 1200^ . Broad-band magnitudes from the 
SDSS-DR7 database we re compared with a grid of the¬ 
oretical galaxy spectra (|Bruzual fc Ghar loti 1200311 span¬ 
ning a large range in star-formation histories. From the 
difference between observed and theoretical magnitudes, 
a likelihood distribution for the mass of each galaxy is 
estimated. The median and the dispersion of this distri¬ 
bution are used for the stellar mass and its error. 

Star Formation Rates (SFRs) are needed for the discus¬ 
sion of the results (e.g.. Figs. [1] and [2). We inf erred them 
from the observed Ha flux as described by iKennicutO 
(|I998ll . The underlying hypothesis is that Ha quanti- 
hes the number of ionizing photons produced by young 
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TABLE 1 

Extremely metal poor galaxies selected in this study 


Name® 

12 -H log(0/H)‘> 

logM*"" 


HII Size® 

AZf 

J020549.3 - 094920.3 

7.58 

± 

0.14 

8.26 

± 

0.48 

26.4 ±1.8 

0.18 

± 

0.03 

Y 

J030331.3 - 010947.1 

7.68 

± 

0.13 

8.33 

± 

0.33 

125 ± 9 

0.72 

± 

0.12 

Y 

J031300.0 -1- 000612.2 

7.58 

± 

0.17 

7.67 

± 

0.31 

119 ± 8 

0.35 

± 

0.12 

Y 

J082555.0 -1- 353231.0 

7.49 

± 

0.11 

6.04 

± 

0.03 

9.6 ±0.7 

0.076 

± 

0.009 

Y 

J094416.6 -1- 541134.4 

7.43 

± 

0.15 

7.05 

± 

0.05 

23.1 ± 1.6 

o.ie 

l± 

: 0.2 

Y 

J113202.4-1-572245.2 

7.59 

± 

0.15 

7.53 

± 

0.48 

22.5 ±1.6 

0.093 

± 

0.021 

Y 

J114506.3 -1- 501802.4 

7.70 

± 

0.10 

6.71 

± 

0.07 

23.6 ± 1.7 

0.31 

± 

0.02 

N 

J115133.3 - 022221.9 

7.32 

± 

0.15 

6.7 

± 

1.0 

12.6 ±0.9 

0.12 

± 

0.01 

Y 

J210455.0 - 003522.0 

7.33 

± 

0.10 

6.19 

± 

0.05 

21.8 ± 1.6 

0.094 

± 

0.021 

Y 

J230210.0 -1- 004939.0 

7.54 

± 

0.12 

8.39 

± 

0.33 

138 ± 10 

0.47 

± 

0.13 

Y 


^ RA and DEC in J2000 coordinates. 

^ Oxygen abundance at the starburst. 

^ Stellar mass of the galaxy from MPA-JHU given in Mq units. 

^ Redshift-based galactocentric distance, in Mpc, from NED (NASA/IPAC Extragalactic 
Database). 

® Seeing-corrected full-width half-maximum size of the starburst, in kpc, as inferred from 
the Hck emission. 

^ Does the galaxy show chemical inhomogeneities? Yes or No. - Note that the only galaxy 
without chemical inhomogeneities is also the object of largest metallicity. 




Distance [arcsec] 


Fig. 2.— Same as Fig. ^ for another example, J0944-I-54. 

(< 20Myr) massive (> 10 Mq) stars, and so, its flux 
scales with the current SFR. 

4. RESULTS 

HII-CHI-mistry provides the abundance of oxygen 
along the major axis of the galaxies. In nine out of 


the ten cases, sharp metallicity drops were found (ex¬ 
amples are given in Figs. [T] and [5J with the only excep¬ 
tion having constant metallicity pointed out in Table [1]). 
Remarkably, the metallicity drops occur coinciding with 
starburst regions (the solid lines in Figs. [Hand[21 bottom 
panels). 

A summary plot with the metallicity of the starbursts 
and the underlying galaxies is given in Fig. [3) The 
abundance of the starburst was determined as the lowest 
abundance, whereas the abundance of the host was com¬ 
puted by averaging all the values in a region of 2.5 arcsec 
around the position with the largest abundance. The 
host galaxies have metallicities higher than the starbursts 
by factors of 3 to 10 where their surface SFRs are smaller 
by factors of 10 to 100. According to the empirically 
established scaling between surface SFR and gas sur- 
face density (th e so-called Kennicutt-Schmidt law; e.g., 
Kennicntilll998t lYang fc Krumhol3l2012t iCormier et alJ 
20141 1. the observed starbursts imply large gas reservoirs 
with surface densities between 10^ and 10^ Mq pc“^ 
(Fig- El upper axis). 

According to simple chemical evolution models the 
starbursts can evolve to become the host only if they 
expel 80% of the gas and metals to the intergalactic 
medium (blue solid line in Fig. [S]). If the evolution pro¬ 
ceeds as a closed-box, and it begins from pure metal- 
free gas, the fraction of mass in metals of the interstellar 
medium Z dep ends only on the gas that remains (e.g.. 
Id mun di[i99g), so that if a region of metallicity Zi and 
gas mass Egasi evolves to reach a metallicity Z 2 of gas 
mass Sgas2, then Zi - Z 2 = Y ln(Egas2/Sgasi) with Y 
the stellar yield, i.e., the mass of new metals eventu¬ 
ally ejected per unit mass locked into stars. The red 
solid line in Fig. [3| shows the evolution expected for a 
region at the red bullet assumin g: a typical Oxygen yield 
(0.004; iMevnet fc Maedei1l2002n . It renders metallicities 
too high to explain the metallicity of the host galaxy by 
close-box evolution of the starbursts. If gas is allowed 
to leak out of the system, the relationship between gas 
mass and metallicity remains formally as for the closed- 
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Fig. 3. — Oxygen abundance of starburst (square symbols) and 
host (asterisks) vs surface SFR for the XMPs with metallicity in¬ 
homogeneities. The axis on top gives the gas surface density for 
a typical gas consumpt ion time scale of one billion years (e.g., 

I Yang Ik Krumholzl 120131 . The lines show the chemical evolution 
of a starburst at the red bullet as it consumes gas and returns 
metals; the closed-box evolution overproduces oxygen (red line) 
whereas leaking out 80 % of the gas explains the difference between 
starburst and host (blue line). 

box model, replacing the stellar yield with an effective 
yield that involves the fraction of mass in stars that 
returns to the interstellar medium i?, as well as the so- 
called mass loading facto r w, that gives the outflow rat e 
in terms of the SFR (e.g.. lSanchez Almeida et al.ll2014al) : 
Yes = Y {l — R)/{l — R-\-w). The blue solid line in Fig. [3 
corresponds to = 0.2 Y, or to w = 3.2 given R = 0.2, 
which implies that for each solar mass of gas locked into 
stars, four solar masses of gas flow out of the system. 

Figure m shows the scatter plot of metallicity vs galaxy 
stellar mass for both the host and the starburst. Our 
objects are outliers of the mass metallicity relationship 
(MZR) observed in the local Universe. The green and 
the blue lines in F ig. ID represent M ZRs for the range of 
low -mass galaxies (|Lee et al.ll200til) and high-mass galax¬ 
ies (jTremonti et al.ll2004[) . respectively, with the shaded 
areas portraying the observed standard deviation at a 
given mass. The high-mass metallicities have been de¬ 
creased by 0.2 dex to make the strong-line method of 
the o riginal work consistent with the DM scale used 
here (jKewlev fc Ellisonl I2008D . The MZR is based on 
the metallicity of the brightest parts of the galaxies 
therefore, XMPs must be represented by the metallic¬ 
ity of their starbursts, which clearly lie below the ob¬ 
served MZR (Fig. |4]). The scatt er at a given mass is 
known to be related to the SFR (jMannncci et al.l [20101 : 
iLara-Lopez et al1l2010t) . so that the higher the SFR the 
lower the metallicity. The pink area in Fig. |4] repre¬ 
sents the spread in metallicity to be expected if the 
SFR changes b y ±10 times the mean SF R given by the 
red solid line (Andrews fc Martliiil I2013D . and this de¬ 
pendence naturally accounts for the large difference in 
metallicity between starburst and host. XMPs are going 
through an active star-formation phase, with starbursts 
of low metallicity. The end of this phase will disclose a 
galaxy having properties close to those of the host, i.e., 
a stellar mass similar to the present XMP bnt with low 
SFR and high metallicity. 


Fig. 4.— Oxygen abundance of the starburst (square symbols) 
and host (asterisks) vs stellar mass of the gal axy. We includ e the 
relationship observed in low-mass (green lin e; ILee et J1I2006I 1 and 
high-mass (blue line; ITremonti et al.ll2004ll local galaxies, the lat¬ 
ter re-scaled as described in main text. The colored green and 
blue regions represent the dispersion at a given mass, therefore, 
XMPs have starbursts that are metal-deficient for their masses. 
The d ispersion of the metal licity vs mass relation correlates with 
SFR IIMannucci et all 1201011 so as to cover the area delimitated 
by the pink region for SF Rs going from ten times la rger to ten 
times smaller than average HAndrews &: Martini|[20T3l . Thus, the 
large difference of SFRs between starburst and host may naturally 
account for their different abundances. 

5. DISCUSSION 

The XMPs studied here have an off-center starburst 
of low metallicity but, otherwise, they are rather normal 
disk-like galaxies with exponential light profiles, mod- 
est rotational velocities, and significant random motion s 
(|Elmegreen et al.l I2012t ISanchez Almeida et al.l 12013D . 
Under these conditions the time-scale for azimuthal gas¬ 
mixing is rather short (see Sect.[T|), therefore, the metal- 
poor gas associated with the starburst must have arrived 
to the galaxy recently. 

The mass surface density at the starburst is unusu- 
ally large for the out er parts of dwarf irregular galaxies 
(iCormier et al.l [20I4H . by a factor of ^ 10. If the ex¬ 
cess mass is from a recently accreted cosmic cloud, then 
it mixed with the ambient gas and diluted the ambient 
metallicity by the same factor, 10:1. This dilution ac¬ 
counts for the metallicity drop in the starbursts. The ac¬ 
cretion time is constrained by the localization of the star 
formation, which is not sheared or spiral-like. The con¬ 
finement of the enhanced star formation to local regions 
several hundred parsecs in diameter (Table [T]) implies 
that the accretion time is shorter than the orbital time. 
For small galaxies like these the orbital time over one 
radian is ^ 100 Myr. Accretion lasting for ~ 100 Myr or 
less corresponds to an impacting cloud with a vertical ex¬ 
tent that is similar to the radius of the galaxy or several 
kpc, considering that the infall speed is likely comparable 
to the rotation speed since both are set by the gravita¬ 
tional potential. From these considerations, we conclude 
the following scenario for XMP galaxies with locally low 
metallicities studied here. Sometime within a period of 
^100 Myr, a cosmic cloud with low metallicity impacted 
the outer region of a regular dwarf galaxy disk, com¬ 
pressing it, diluting the original metals, and triggering 
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excess star formation, all by factors of 10. This picture 
is consistent with the current ideas on h ow dwarf galaxies 
grow through gas accre t ion ('Sect.fTlandlFinlator fc Dav^ 
l2f)M IDekel enil[2?M [Scha.ve et ahlMToir even though 
the interaction between gas streams and forming disks re¬ 
mains to be proper ly understood (jCeverino et al.l [20101 : 
iVerbeke et al.l[20Tll . 

Explanations alternative to the cosmic origin of the gas 
cannot be completely ruled out, although they are dis¬ 
favored by observations. Major mergers are not obvious 
in the optical images, and can be further discarded argu¬ 
ing that XMPs reside in low density environments, with 
their large HI envelops showing only moderate distor¬ 
tions (|Filho et al.l 1201^ . Minor gas-rich mergers could 
explain localized metallicity drops, provided that the gas 
in the accreted galaxy has a metallicity lower than the 
metallicity at the drop. A minor-merger scenario dif¬ 
fers from a gas accretion scenario only if the merging 
galaxy has enough stars and dark matter. The addi¬ 
tional mass would generally be much larger than the 
merging gas mass alone, and the distortion on the host 
would be larger as well. Merging stars should also pen¬ 
etrate the host disk, unlike the merging gas, and these 
stars would orbit differently than the host stars, possibly 
producing star streams in the host halo. Since neither 
large HI distortions nor stellar streams are evident in 
XMPs (see, e.g., Figs.[T] and[2]), the hypothetical merging 
galaxy must have little dark matter and stars and, conse¬ 
quently, it could be likened to a gas stream. On the other 
hand, model cosmic gas streams are often c lumpy, with 
some stars forming in the densest cores fe.g.. lDekel et al.l 
[200l), in which case the accretion of a clumpy stream is 
equivalent to a merger with a small very gas-rich dark- 
matter-poor satellite. The difficulty to distinguish a pure 
gas accretion event from a gas-rich minor merger can be 
linked to the problem of finding an operational defini- 
tion of galaxy at the low-mass end of galaxy mass es (see 
iForbes fc Kroupal 1201 H iWillman fc StradCT 1201211 . Yet 
another possibility for feeding XMPs could be gas accre- 
tion from the galaxy halo i nduced by galac tic fountains 
(jFraternali fc Binnevll200§ : lFraternaljl2ni4l l . The metal- 
rich gas ejected from the stellar disk acts as a catalyst 
when it becomes mixed with the hot metal-poor halo gas. 
The mixing catastrophically decreases the halo cooling¬ 
time, and mass from the halo is dragged along when the 
fountain gas falls back onto the disk. Several observa¬ 
tions disfavor this possibility, though. The mass dragged 
along is only a fraction of the original mass, and so, the 
gas falling in is not very metal-poor. In addition, the 
accreted material is expected to spread out over the disk 


rather than being concentrated, and the gaseous haloes 
in dwarfs are not massive, suggesting the mechanism to 
be inefficient in our particular case. 

A unidirectional propagation of star-formation may 
produce a head-tail morphology, with decreasing surface 
brightness and increasing age and metallicity from the 
position of t he current starburst. T his mechanism was 
proposed by iPanaderos et al.l (1199811 to explain the color 
gradient in the XMP galaxy SBS 0335-052E, and may 
work for other XMPs as well. The required unidirec¬ 
tional propagation would naturally occur if the under¬ 
lying gaseous structure is already very prolate. Even if 
uncommo n, cigar-like galax ies are theoretically conceiv¬ 
able fe.g.. lBaek et al.l 1200311. and hav e been observed at 
high redshift (e.g., lYuma et al.ll201lll . 

The gas-cloud impact scenario explains the lower limit 
for the metallicities of the XMPs, which is around l/50t h 
the solar metallicity fe.g.. lSanchez Almeida et al.ll2ni4^ . 
The absence of objects with lowe r metallicity has chal¬ 
lenge d interpretation for decades (iKunth fc Lebouteilled 
[2011D. However, the threshold is naturally set in this 
context by the metallicity of the cosmic web gas at 
the present time, wh ich is predicted to be around the 
observed limit (e.g., iFumagalli et al.l [201111 . The sce- 
nario also explains why XMPs tend to be co met ary 
(|PaDaderos et al.ll2008t iMorales-Luis et al.l [201111 . The 
impacting gas clouds have the largest effect on the out¬ 
skirts of galaxies where the ambient pressure and column 
density are low. For most orientations of the galaxy, this 
gives a cometary shape (like Figs. [Hand[2]). 

If our scenario is correct, XMPs should not be regarded 
as pristine galaxies but as objects that are still forming 
in the nearby Universe. Their gas traces the intergalactic 
medium in their immediate surrounding, providing a new 
method to characterize the cosmic web. 
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